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Abstract 



■ ■ The paper contains mathematical justification of basic facts concerning the 

Q,' Brownian motor theory. The homogenization theorems are proved for the Brownian 

motion in periodic tubes with a constant drift. The study is based on an apphca- 
tion of the Bloch decomposition. The effective drift and effective diffusivity are 
expressed in terms of the principal eigenvalue of the Bloch spectral problem on the 
cell of periodicity as well as in terms of the harmonic coordinate and the density of 
the invariant measure. We apply the formulas for the effective parameters to study 
the motion in periodic tubes with nearly separated dead zones. 
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1 Introduction 



K^ I The paper is devoted to mathematical theory of Brownian (molecular) motors. The 

concept of a Brownian motor has fundamental applications in the study of transport 
processes in living cells and (in a slightly different form) in the porous media theory. There 
are thousands of publications in the area of Brownian motors in the applied literature. For 
example, the review by P. Riemann [6] contains 729 references. Most of these publications 
are in physics or biology journals and are not mathematically rigorous. Some of them are 
based on numerical computations. 

Consider a set of particles with an electrical charge performing Brownian motion in a 
tube Q with periodic (or stationary random) cross section (see Fig. 1). 
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Figure 1: A periodic tube Q with "fingers". One can expect that Vcg < \V^q\. 

We will assume that the axis of Q is directed along the xi-axis. Let's apply a constant 
external electric field E along the axis of Q. Then the motion of the particles consists of 
the diffusion and the Stocks drift, and the corresponding generator has the form 

ox I 

complemented by the Neumann boundary condition on dQ (we assume the normal reflec- 
tion at the boundary). 

One can expect that the displacement of the particles on the large time scale may 
be approximated (due to homogenization) by a one-dimensional diffusion process Xi{t) 
along the xi-axis with an effective drift Vcg and an effective diffusivity a^, which depend 
on V and the geometry of the tube Q. If we reverse the direction of the external field 
from E to —E, then the corresponding effective drift V^ and the corresponding effective 
diffusivity (cr^)^ will be, generally, different from V^eff and a^ when Q is not symmetric 
with respect to the reflection xi — )■ —Xi. For example, one can expect that V^s < |V^il 
for Q shown in Fig. 1. The difference between the effective parameters can be significant. 
Then by changing the direction of the exterior electric field E periodically in time, one 
could construct a constant drift in, say, the negative direction of xi and create a device 
(motor) producing energy from Brownian motion. 

The idea of molecular motors goes back to M. Smoluchowski, R. Feynman, L. Brillouin. 
Starting from the 1990-s, it became a hot topic in chemical physics, molecular biology, 
and thermodynamics. The following natural problem must be solved by mathematicians: 

1. The homogenization procedure (reduction to a one- dimensional problem) must be 
justified. 

2. Expressions for effective parameters V^g , a^ in terms of an appropriate PDE or a 
spectral problem on the period of Q (in the case of periodic tubes) must be found. 

3. For some natural geometries of Q that include a small parameter, asymptotic 
expressions for effective parameters need to be obtained. 

We solve the first two problems here using an analytic approach that justifies the 
homogenization procedure and allows us to express the effective parameters in terms of 



the principal eigenvalue of a spectral problem on one period of Q. We will also provide a 
couple of simple consequences of the obtained formulas. We show that dVcs/dV = cr^ > 
when ^ = 0. Thus, since Vcs is analytic in V, it is a strictly monotone function of V 
when 1^1 is small enough. Simple asymptotic formulas for the effective parameters will 
be justified for periodic tubes with nearly separated dead zones. Note that the first two 
problems listed above also can be solved with probabilistic techniques similar to those 
used for homogenization of periodic operators in M"', as we'll discuss in a forthcoming 
paper. 

Our approach is based on the Bloch decomposition. Recently several papers appeared 
(see [2], [3] and references there), where the Bloch decomposition was used to study ho- 
mogenization problems in periodic media. It was shown that this approach has many 
advantages (when it is applicable). So far, this approach was applied to self-adjoint ellip- 
tic equations in the whole space or in a domain with a finite boundary. The symmetry 
and the existence of a bounded inverse operator were essential there. 

We consider a parabolic problem. It is non-self-adjoint due to the drift, which can not 
be neglected since it is essential for applications. Besides, the homogenization is applied 
only with respect to one variable. All these features (the non-symmetry is the main 
difficulty) make the problem under consideration essentially different from the applications 
of the Bloch decomposition in the homogenization mentioned above. 

We also would like to note the presence of the boundary integral term in the formula 
for the effective drift (see Theorem 12.51) . In probabilistic terms, it arises from the re- 
normalized time spent by the Brownian motion with the drift on the boundary (local 
time). Our next paper based on a probabilistic approach to the problem will contain the 
detailed analysis of the asymptotic behavior of the effective parameters with respect to 
the parameters describing the geometry of the domain. 

The plan of the paper is as follows. Lemma fI72\ describes the properties of the principal 
eigenvalue Aq = Ao(6') of the Bloch spectral problem on the cell of periodicity of Q. The 
main results on homogenization are obtained in Theorems 12.11 and 12.31 In particular, it 
is shown there that the effective drift V^s and the effective diffusivity a"^ are given by the 
coefficients of the Taylor expansion of the eigenvalue Ao(^). To be more exact, Xo{9) = 
iVeffd — a'^6'^ + 0\6\^, ^ — 7- 0. Expressions for the effective parameters through the harmonic 
coordinate and the invariant measure are given in Theorem 12. 5[ The latter formulas are 
applied to a particular class of domains in the last section. Periodic tubes with nearly 
separated dead zones are considered there. We rigorously justified the asymptotic formula 
for the effective diffusivity, which was found earlier for domains with somewhat simpler 
geometry in pQ, [1]. 

2 Description of the model and main results. 



Consider a tube Q C M*^, d >2, periodic in xi with period 1, with a smooth boundary dfi 



see Fig. 1). Denote by S^-^^ the cross-section of Q by the plane M.'^ ^ orthogonal to the 



Xi-axis at the point Xi We assume that Sx^ is bounded (and periodic with respect to xi). 

In the simplest case, the boundary of the cross-section is a function of the angular 
variables. For instance, Sx^ is the ball of the radius R (xi) if dVt is the surface of revolution. 
However, in general it may have a very complicated form. For example, Vt in Fig.l contains 
"fingers" ("dead ends" in the terminology of [7], [8]), and Sxi is not connected when 
a < xi <h. 

Consider the cell of periodicity in Vt defined as Vt' = f^HlO < a^i < !}• It is assumed 
that the compact manifold obtained from Q' by gluing 5*0 to 5*1 is connected. For trans- 
parency, one can assume that Q' itself is connected. For example, Q' in Fig.l is connected 
if the origin is slightly to the left of a. It is not connected if the origin is the midpoint be- 
tween a and b, however, it still forms a connected manifold after gluing 5*0 to 5*1. One also 
can define a connected cell of periodicity Q' in Fig.l using planes M.'^^^ through (a -|- b)/2 
and (a + 6)/2 + 1 by cutting only the central narrow part of Q, but not cutting the fingers. 
In the latter case, we use the notation 5*0, Si not for the whole cross-sections of Q, but 
for the parts of the boundary of Q' that belong to the corresponding M*^"^ planes. 

Denote the exterior normal to dQ by n. Let x = {s, z), where s = xi, z = {x2, ...,Xd)- 
Consider the following parabolic problem in the tube Q: 

OZL OZL OZL 

— = Au + V—, X eQ; ^|af7 = 0; u{0,x) = if{es,z),xeQ, (1) 



where ip G Li^com(IR'^) (the space of integrable functions with compact supports). 

The main results of the first part of the paper are stated in the following two theorems. 
Theorem 12. II specifies the asymptotic behavior, as e J, 0, of the solution u of the problem 
([1]) in terms of the solution w of the one- dimensional problem 

-KT = a -— + Vcs-^, s EK; w\t=o = wo{s,es). (2) 

Here cr^ is the effective diffusivity, Vcs is the effective drift, and the initial function Wq{s, es) 
Remove the arguments? is an average oi(f in 2;- variable and will be defined later in (|5]). 
This theorem provides also the main term of asymptotics of u in terms of the solution W of 
a one-dimensional, e-independent problem if ip is smooth. The second theorem describes 
the effective parameters a^ and Vcs in terms of the problem on the cell Q'. 
We will use notation 0{e) for functions / such that |/| < Ce, e J, 0. 

Theorem 2.1. 1) Let ip G Li^comO^''') ■ Then for each to > 0, 

u = w + 0{e), t > —, 

where w is the solution of ^ and Wq is defined in ^. 

2) Let if and its derivative ips he continuous hounded functions with compact supports. 
Let t = ^, s = ''" "^'^ and let u(t, q, z) = u{t, s, z) he the function u written in the new 
variables. Then for each tq > 0, 

u = w + 0{e), T > To, 



where w is the solution of the following e-independent problem 



dw r,d w 



aTTT, /■! 



a 



2' 



dr dq"^ 



^ G M; w\t=Q = ^{s)= / wo(s,?)rfs 







Remarks. 1. An explicit form of Tp is given below in (|5]). 

2. The proof of the theorem allows one to obtain explicit estimates of the reminder 
terms in the formulas above through the corresponding norms of ^p . 

We need to state one important lemma before we can formulate the second theorem. 
Denote by S^ the lateral part of the boundary dO! of ^2', i.e., 5' = 9r2'\[S'o IJ Si]. Consider 
the following eigenvalue problem on the cell Vt': 

(11] (jlj 

/^v + V— = Xv, X e n'; —y = 0; ^;|5, = e^^Uo, <\sr = e^'iQso- (3) 

This is an elliptic problem in a bounded domain, and for each 6 the spectrum of this 
problem consists of a discrete set of eigenvalues A = \j{0), j > 0, of finite multiplicity. 

Lemma 2.2. Let 9 G [— 7r,7r]. There exists a simple eigenvalue A = Ao(6') of problem ^ 
in a neighborhood \9\ < 5 of the origin ^ = and real constants Kfj, a such that 

1) Ao = iV^^e-a'^e'^ + 0{e^), ^ -> O, where a^ > 0, 

2) ReAo(^) < 0, < 1^1 < 5, 

5; ReA,(e) < -71, \e\<5, j>0, 71= 71(5) >0, 

4) Re\j{e) < -72, 5<\e\< vr, j > 0, 72 = l2{5) > 0. 

Theorem 2.3. The effective diffusivity a^ and effective drift Vcs defined in ^ coincide 
with the constants introduced in Lemma \2.2[ 

Let us describe the averaged initial function if. We need the problem adjoint to 
problem ([3]). This problem has the form: 

Av-V— = Xv, X G n'; (— - Vni)v\si = 0; v\s, = e'%\so, v'sls, = e''vi\s,, (4) 

where ni = ni{x) is the first component of the normal vector n. Then A = Ao(6') is an 
eigenvalue of problem (j4]). Let v = iPq{9,x) be the corresponding eigenfunction, and let 
tt{x) = z^q(0,x) be the eigenfunction of (j4]) when ^ = (and Aq = Ao(0) = 0). Then Wq,^ 
are defined as follows: 

Wo{s,es)= / 7r{s, z)ip{es, z)dz, Ip = / n{s, z)(p{(^, z)dzds, (5) 

Jss Jo Jss 

where Sa is the cross-section of Q by the plane through s = a. 



The following lemma will be needed in order to prove the theorems above. Consider 
the non- homogeneous problem ([3]) in the Sobolev space H^{S): 

(A + V- - A).; = /, X G n'- —y = 0; v\s, = e%\s,, <U. = e'^slso- (6) 
This is an elliptic problem, and the resolvent 

is a meromorphic in A operator with poles at a discrete set of eigenvalues A = \j{9). 
Denote by Qa the sector in the complex A-plane that does not contain the negative semi- 
axis and is defined by inequalities — vr + a < argA < vr — a. 

Lemma 2.4. For every a > there exist constants R = R{a) and C = C{a) such that 
the region Qa,B. = Qaf]{\M > -^} does not contain eigenvalues \j{9), and the following 
estimate is valid for the solution v = R^^f of the problem ^: 

\M\ML,in')<C{a)\\f\\L,^n'), >^ e Q^,r. (7) 

Proof. This lemma can be proved by referencing the standard a priori estimates for 
parameter-elliptic problems [5j. We will provide an independent proof. We multiply 
equation ([6]) by w and integrate over Q'. This leads to 

- / {\Vv\^ + X\v\'^)dx + / Vv'^vdx= / fvdx. (8) 

Jn' Jw Jn' 

The terms on the left in ([8]) can be estimated as follows 

I / {\Vv\^ + \\v\^)dx\ > ci(«)(||Vt;||i^(f,,) + |A|||t;||i^(f,,)), A G g„, 
Jn' 



and 

^^' /||2 *' IIII2 

V VgVUX\ ^ ||(' 



ci(a)„ /„2 , V^ 



Vv'^vdx\ < ^— ||<||l2(o') + ^^;^ll^ll^2(f^')- 



We put R = , y^-)i2 - Then the absolute value of the left-hand side in ([8]) is estimated from 
below by ^2 (Il^^lli2(^') ~'' I'^lll^'lli2(n')) '^^en A G Qa,R, and (^ implies that 

^(l|Vt;||i,(f,,) + |A|||^;||i,(^,)) < \\fKin')\\v\\L,in'), A G Q„,r. (9) 

We omit the first term on the left and obtain (jTj). D 

Proof of Lemma 12.21 Denote hy Hq = A + V^ the operator in L2(fi') defined on the 
functions from the Sobolev space H'^{Q') satisfying the boundary conditions (E]). One can 



define the parabolic semi-group e*^". Its integral kernel is the fundamental solution of 
the corresponding parabolic boundary value problem in Q'. When 6 = 0, this kernel is 
real and positive. Thus Perron-Frobenius theorem is applicable to the operator e*^", i.e., 
the operator e^^° has a unique maximal eigenvalue /i such that /i > 0, // is simple, the 
corresponding eigenfunction is positive, all the other eigenvalues fij (perhaps complex) are 
strictly less than /i in absolute value {\fij\ < /i), and the operator does not have strictly 
positive eigenfunctions with eigenvalues different from /i. 

We note that Uq = 1 is an eigenfunction of i^o with the eigenvalue Aq = 0. Thus vq is an 
eigenfunction of e*^° with the eigenvalue fi = 1. Since wq > 0, from the Perron-Frobenius 
theorem it follows that /i = 1 is the maximal eigenvalue of e*^°. This implies that 
ReXj < Aq = for all the eigenvalues \j 7^ of the operator Hq. Since the eigenvalues 
Xj form a discrete set, from Lemma 12.41 it follows that there exists a 7 > such that 
ReXj < -7 for all A^ 7^ 0. 

Since Hq depends analytically on 6, the location of the eigenvalues Xj when ^ = 
and Lemma 12.41 imply that there exists 6 > such that the operator Hg has the following 
structure of the spectrum when 16*1 < 6: the operator has a simple, analytic in 6 eigenvalue 
A = Xq{9), Ao(0) = 0, and ReAj(6') < —71 = —7/2 for all other eigenvalues Xj{9) of Hq. 
The statement 3 of Lemma 12.21 is proved. Consider now the eigenvalue A = Ao(^) for 
purely imaginary 6 = iz, z > 0. li v belongs to the domain of the operator Hi^, then v 
also belongs to the domain of the operator Hi^. From here it follows that both A = XqIiz) 
and A = Xo{iz) are eigenvalues of Hiz. Since, A = Ao(^) is the unique eigenvalue in a 
neighborhood of the point A = 0, it follows that Xo{iz) is real. This implies statement 1 of 
Lemma [2.21 except the inequality a^ > 0. The latter inequality will be proved below (in 
Theorem 12. 5p by a direct calculation of a"^ (see formula (126!) ). Thus it remains to justify 
statements 2 and 4. This will be done if we prove that for every 61 > 0, there exists 
7 = 7(^1) > such that ReXj (9) < —7 for all the eigenvalues of the operator Hq when 9 
is real. Si < \9\ < n. 

Let us prove the latter estimate for the eigenvalues Xj{9), 61 < \9\ < ir. We fix 
an arbitrary rational 9 = 9' = ^^ in the set 5i < \9\ < vr. Consider the domain 
r2 = fi f]{0 < Xi < n}, which consists of n elementary cells of periodicity D,^^"^ = f2 f]{j < 
Xi<j + 1}, < j < n — 1 (n(°) coincides with the previously introduced cell fi'). The 
lateral side of the boundary of the domain Q will be denoted by S'-, and the parts of the 
boundary located in the planes through the points Xi = and Xi = n will be denoted by 
5*0 and S'„, respectively. Let Hg = A -|- V^ be the operator in L2(n) that corresponds 
to the following analogue of the problem (j3]): 

^v + V—= av, xeQ; — |g, = 0; v\s„ = e^^l^o, <|5„ = e^"'<|5o- (10) 

We denote its eigenvalues by a = aj, and we keep the notation A = Xj for the eigenvalues 
of this operator when n = 1 (and Q = Q'). 

When 9 = 9', the conditions on Sq, Sn become the periodicity condition. Thus the 
spectrum of Hq? has the same structure as the spectrum of Hq, i.e., ao = is an eigenvalue 



with the eigenf unction vq = 1, and 

Rea^ < -7 < (11) 

for all other eigenvalues aj of Hg? . 

We compare the set {Aj} of the eigenvalues of the operator Hg' and the set {aj} of 
the eigenvalues of the operator Hg, . The following inclusion holds: {Aj} C {aj}. Indeed, 
if V is an eigenfunction of the operator Hgr, then one can construct the corresponding 
eigenfunction of Hg, with the same eigenvalue by defining it as e^^^v in each elementary 
cell il^^^ Cf2, 0<j<n— 1. However, these two sets of eigenvalues do not coincide. 
In particular, A = is not an eigenvalue of Hgi (while a = is an eigenvalue of Hq,). 
Indeed, from the simplicity of the eigenvalue a = it follows that A = could be an 
eigenvalue of Hq/ only if vq = 1 is its eigenfunction, but Vq does not satisfy the boundary 
conditions ([3])- Thus, fITT]) implies that ReAj < —7 < for the set of eigenvalues of the 
operator Hg/. Then the same estimate with 7/2 instead of 7 is valid for the eigenvalues of 
Hg when ^ is in a small enough neighborhood of 9' = ^^. One can find a finite covering 
of the set {6 : 61 < \6\ < ir} by some of these neighborhoods. Thus, the desirable estimate 
of eigenvalues Xj{6) is valid for all 6 of the set {^ : ^i < \6\ < tt}. 

D 

Proof of Theorems 12.11 and 12.31 The solution u of problem ([1]) can be found using 
the Laplace transform in t: 

fa+ioo 



u = — :/ v{X,x)e^^dX, a > 0, 



where v G Iv2(fi) is the solution of the corresponding stationary problem: 

dv dv 

^v + V—-\v = -v9, xen- —\an = 0. (12) 

OS on 

In order to solve 012p . we apply the Bloch transform in the variable s = Xi: 

00 
w{x) -^ w{6,x) = '^w{s -n,z)e'"'^, 9 E {-it, it), x e ^. 



This map is unitary (up to the factor l/A/27r) from L2(fi) to L2([— 7r,7r]) x L2(fi'), and 
the inverse transform is given by 



1 r - 

w{x) = — / w{9,x)d9, X E Q. 



Note that function e '^^^w{9,x) is periodic in s. Thus the knowledge of w on Q' allows 
one to recover w on the whole tube fi. 

8 



The Bloch transform reduces problem (fT2l) to a problem for v on the cell of periodicity 
Q'. The latter problem has form ([6]) with / = —if, i.e., v = —R^^(p. Thus 



u(t,x) 



Air'^i 



a+«00 PTT 



yS ,^„Ai 



Rl^pe^'dddX, xen. 



(13) 



Function Rf^ip here is defined in Q', but it is extended to the whole tube Q in such a way 
that e~'^^'^R\(p is periodic in s. 

Now that formula flT3|l for u is established, we are going to study the asymptotic 
behavior of m as e — )■ 0. We split u as u = ui + U2, where the terms ui,U2 are given 
by flT3l) with integration in 6 over sets |^| < 5 and 6 < \6\ < vr, respectively. We choose 
6 small enough, so that Lemma 12.21 holds, and then make it even smaller (if needed) to 
guarantee that — ^ <ReAo(^) < 0, \9\ < 5. 

Im X 



r. 



ReX 



Figure 2: Contour F = Fi + r2 + F3, where 7 coincides with 71 or 72 defined in Lemma 
12. 2[ and F2 is long enough, so that the estimate ([7]) is valid on Fi and F3. 

Let us show that M2 does not contribute to the main term of asymptotics of u. Lemmas 
12.41 and 12.21 allow us to rewrite U2 in the form 

U2(t,x) = ^ Rl(pe^*ded\ x G fi, (14) 

47r^« Jt J5<\e\<n 

where F is the contour shown in Fig. 2 with 7 = 72. From Lemmas 12.41 and 12.21 it follows 
that 



Cl 



WRI^Wmw) < --j^ML.iw), A G F, 5<\9\< 



TT. 



Thus 



s<\e\<w 



\\Rl^\\L,in')d9 < - 



Cl 



+ A 



|¥'||L2(n')xL2(|e|<7r) 



C2 



■\\^{es,z)\\L2^n) 



< 



C2 



1+|A| 

(1 + |A|)V^"^^'^)II^^(^'')' ^^^- (^^) 

One can replace fi' here by a shifted cell of periodicity fi*^-'^ defined by inequalities j < 
Xi < J ' + 1, since the function e~'^^^^R\(^ is periodic in xi. We put (IT5ll with fi' replaced 
by fi^-'^ into f lT^ and estimate the integral over F. This gives 

Ik2(t,x)||i2(j^0)) < — =e-^*||v5(x)||i2(Rd), 



where c does not depend on j. 

Now one can obtain a uniform estimate of U2 by using the following inequality for 
the Green function G{t, x, y) of the problem ([T]): |G(1, x,y)\ < ae"^*-^^"^^-* , which implies 
that 

f .00 

\u2{t,x)\< / \G{l,x,y)u2{t-l,y)\dy<^e-^' ^ e-''(^-^-)'||^(x)|U,(i,.) 

■^^ ^^ i=-oo 

C 

< ^e"'^*||v9(x)||L2(]8d), X eQ, t>l. 

In particular, U2 = 0{e) when t > to/e^. Hence U2 does not contribute to the main term 
of the asymptotics of u. 

In order to study ui, we also shift the contour of integration (in A) to F. Now we take 
7 = 7i, and in this case the simple pole of i?^ at A = Ao(^) must be taken into account. The 
residue at this pole of the integral kernel of the operator i?^ is equal to —il'{9, x)ip*{0, y), 
where ip is the eigenfunction of the problem ([3]) and ip* is the eigenfunction of the problem 
(j4]) (the residue is a negative operator). Thus 

^'^J\e\<sJn' 47r^WrJ|e|<5 

(16) 
where the functions ipiip* are extended to ^2 using the Bloch periodicity condition, i.e., 

^{e,x) = e'''ij,{9,x), r{0,x) = e^'^roiO^x), (17) 

where the functions V'o; V^o ^^^ periodic with respect to s. The second term in fITB]) can be 
estimated similarly to U2. Thus 

u{t, x) = — I I tp{e, x)if{e, y)e^'''-^^^i^dyde + w{t, x), 
27r J\e\<s Jn' 

where x G fl, t > 1, and 

\w\ < -^e-^*||<^(a:)|U,(M.) = 0{e) ii t > %, ^ e Li,,^^. (18) 



Let y = {s',z'). Then 



^ 



J2v{e{s'-n),z')e^^'. 
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Hence (with (IT71) taken into account), we have 

^(^' ^) = ^y. I I e''(^^^'+"Vo(^, xWo{0, ?/)e^°WV(^(5' - ^), ^')dyd6 + 0(e) 
= ^E /" / e^'(-^'Vo(^,a;)V^(^,2/)e^"WV(«', ^0^2/^^ + 0(e) 

3*^^'^''Vo(^, a;)^(^, 2/)e^«(^)V(£s', /)ciyci^ + 0(e).(19) 



27r ./ifl|<5 jfi 



We spht the double integral / above in two parts I = Ii + 12, where the integration in 
6 in the part Ji extends only over the segment [^"^tl < 1, t > 6^^, and I2 = I — h- From 
Lemma [2.21 it follows that ReAo(^) < —^6^, \6\ < 6, ii 6 is small enough. Since we can 
take 6 as small as we please, and the functions '^/'oj'^o ^'^^ bounded, it follows that 



h<C I I e~'r<^''\^^es',z')\dyd9 = }}Z}lIltl^o{e-'T^''') = 0{e), t > ^. 



We replace Xo{9) in the integral Ji by its quadratic approximation found in Lemma [2^ 

If only the remainder term is taken into account, then the corresponding integral can be 
estimated by 

C f [ \9\He-'^''\{es',z')dyd9 < ^y{y)\\^^^^,^ T \9\He--''''d9 

>/|0|<t-i/3 Jn s j„oo 

From here and (IT9l) it follows that 



■My)\\L^(K-') = 0{e), ift>^2' 



^(t^^) = ^l I e''^^^^'^MO,x)ro{9,y)e^''''^''^-'''Mes',z')dyd9+0{e), t > ^. 

'\e\<t-^/'-^ Jn 



^TT ./|fl|^+-l/3 ./o £ 



We put ^ = in the arguments of functions ipQ and -j/^q. Since these functions are 
periodic, ipQ = 1 and ipliQ.y) = Tx{y) is real, we obtain that ipo{9,x)ipQ{9,y) = 7i{y) + 
0(|6'|). We put this relation into the formula above and note that the integral with the 
term 0(|6'|) does not exceed 



C [ [ \9\e-'^"'\{es',z')dyd9 < ^y{y)\\^^^^,^ f 

J|e|<i-i/3 Jn ^ J-c 



r"00 

\9\e-''"'"'d9 



^"--' 0{e\ ift>5. 
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Hence if t > tq/s^, then 

u{t, x) = ^ f [ e'''^'-''^7r{y)e'^'^^«'-'''''^'ip{es', z')dyde + 0(e), y = (s', z'). 



27r j|e|<t-i/3 Jn 

We can replace here the integration in 6 over the interval 16*1 < t~^''^ by the integra- 
tion over the whole line since the difference between the corresponding integrals decays 
exponentially as i — ?■ oo. Thus 

u{t, x) = — r f e''^'-''\{y)e^^^^'-''''''^'^{es\ z')dyde + 0(e), t > To/e\ 

The integral w = w(t,s) above does not depend on z. By simple differentiation one 
can check that the integral satisfies equation ([2]). Function w{0, s) is the Fourier transform 
in s' followed by its inverse (in 6), i.e., 

w{0, s) = — / / e'^(^-^')7r(s', z')ip{es', z')dz'ds'de = / 7r(s, z')^{es, z')dz', 

(20) 
where Sa is the cross-section of fi by the plane through s = a. Thus w is the solution of 
problem ([2]) with wq given in ([5]), and the first statement of Theorem 12.11 is proved. 

In order to prove the second statement of Theorem 12.11 we will show that w — w = 
0(e), r > tq. We denote function w in new coordinates by w, i.e., w{t, q) = w{^, ''~^°^'^ ). 
Obviously, w satisfies the same equation as the equation for W, and w{0,q) = u;o(-,^). 
The solution of the initial problem for w is the convolution of the smooth (when r > tq) 

kernel K = J- — e^sr and wo(-,?). The second factor is periodic in the first argument 
(see ([5])). Thus, the convolution differs by 0(e) from the same convolution when the 
second factor is replaced by its average with respect to the first argument. This second 
convolution is w. This completes the proof of the second statement of the Theorem 12.11 

It remains to note that in the proof of Theorem 12.11 it was shown that the effective 
drift and the effective diffusivity were defined by the coefficients of the Taylor expansion 
of the eigenvalue Ao(^). Thus, Theorem 12.31 was also established. 

D 

In conclusion of this section we will provide some formulas that can be useful for 
practical evaluation of the effective drift Vcff and the effective diffusivity o"^. 

Let 

v = \^iQvx{x)-Q'^V2{x)^0{B^), e^O, (21) 

be the Taylor expansion at zero of the principal eigenfunction of the problem ([3]) with the 
eigenvalue A = Ao(6'). We plug fl^ and the expansion for Ao(6') from Lemma fl?2\ into (^ 
and take into account that conditions on 5*0, 5*1 in (^ can be rewritten as the periodicity 
in s = Xi of the function e^^^v. Then we obtain the following problems for the coefficients 

At;i + V^ = Ves, X e n'; -^y = 0; vi = s + ^i(x), (22) 

as on 
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ill ) /ii ) Q 

^V2 + V—^ = V,sVi + a\ xen'; -t^\s' = 0; V2 = - + sM^) + M^), (23) 
OS on 2 

where the functions 'ipi,ip2 are periodic. Function vi is called the harmonic coordinate. 
Usually the harmonic coordinate satisfies a homogeneous equation, but this equation 
(see ( l22l) ) becomes inhomogeneous in the presence of a drift in the problem. Denote by 
IT = 7r(x) the principal eigenfunction with eigenvalue A = for the adjoint problem (jll) 
with 6 = 0, i.e., 

An-V— = 0, X eCl'; {- Vni)n\si = 0; n\s, = n\s,„ n'^\si = rr'Jso- (24) 

This is a positive function due to the Perron-Frobenius theorem (see details in the proof 
of Lemma [23]) • We normalize vr by the condition J^, ndx = 1. This function is the density 
of the invariant measure. Note that 7t{x) = l/\^'\ ii V = 0. 

Theorem 2.5. The effective drift V^s can be found from either of the following three 
formulas: 

Ves= f{ViT-7i',)dS = V- f ■K'^dx = V- f nmdS, (25) 

Js Jn' Jds' 

where S is a cross-section of the domain Q by an arbitrary hyperplane s = const., n^ is 

the first component of the outward normal vector n, and 5' = 9r2'\[S'o IJ Si] is the lateral 

part of the boundary of Q! . The effective diffusivity is given by 

a^ = / \Vvi\^^dx. (26) 

Remark. If l^ = 0, then formula ( l26l) is equivalent to the following one 

The latter formula shows that the effective diffusivity is always smaller than in the free 
space if Vt is not a cylinder (i.e., if ifji is not identical zero). 

Proof. The remark will be justified at the end of the proof of the theorem. Let L = 
A + V-K^ be the differential expression in the left-hand side of fl22l) . and let L* = A — V^ 
be the conjugate expression (see (^^). Since L*tt = and Vi satisfies ( 12^ . we have 

{Lvi)Txdx — / viL*Txdx = Kfr / T^dx = Vcs- 
Jq' Jn' 

The left-hand side here can be rewritten using the Green's formula and the divergence 
theorem applied to the vector field F = (t>i7r, 0, 0). This leads to 

Kff / irdx = / [{Avi)7v -ViAn + V-^^^^^^]dx= / {^ir - Vi^ + VviTrni)dS. 
Jn' Jaw ^-^ Jdn' ^'^ "^^^ 
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The second integrand vanishes on S'' since vr satisfies the boundary condition ( 124|) on S'' 
and fi = on SK Thus 

on 

Veff f 7rdx= f (^TT - t;i^ + Vv^^Tn^)dS. (27) 

We substitute here Vi = s + ipi. The integral with ipi is zero due to the periodicity of the 
functions ipi and tt. Thus flTTI) holds with t^i replaced by s. If we also take into account 
that 

/ ^^dS = 0, 

then we obtain the first equality ( 125|) with S" = 5*1. Then this equality holds with any 
5" since V^s is invariant with respect to the shift s — )■ s + a. Let us provide another 
way to show that the middle term in ( 125|) does not depend on the choice of S. Let 
fia,b = fin{«< s < b}- Then 



= / L*7rdx = / (tt^ — 1^7rni)(iS' = / (tt^ — V^7rni)(iS' 

= I {<- VTx)dS - / « - VTl)dS. 

JSb JSa 

In order to prove the second equality i^^, we write the first one with S = Sa and 
integrate with respect to a over the interval (0, 1). The last equality (125]) follows from the 
divergence theorem. 

The proof of (126|1 is similar. We have 



/ (VcsVi + a^)TTdx = / {Lv2)TTdx — / V2L*7rdx = / (t^— ^ ~ ^'27^ — I- V^f27™i)(iS'. 

Jn' Jn' Jn' JsoUS^ dn dn 

2 

We plug here f 2 = y + sipi + i/j2- The integral with 1IJ2 vanishes due to the periodicity of 

2 

ip2 and TT. Hence f2 can be replaced by ^ + sipi. We reduce the integral over So to an 
integral over Si using the substitution s — )• s — 1. Taking into account the periodicity of 
ipi and 71 we obtain 



2\_7 / r/1 , ^''/'1n_ /I , / N^^ , TrA 



j {VeEVi + ^')vrrfx = / [(1 + ^)7r - (| + ^1)^ + V{^ + V'i)7r]d5 



/■ .dvi d-K ^^ ^ ,^ 1^^ 



The last relation is the consequence of the first equality fl25|) and the equality i^i = s + ipi. 
Hence, 



/ [^TX -vi— + VviTi]dS -V^s j ViTldx--Vc{ 



a = [^^TT — 'L>i^:^ — \- VviTr]dS — Vcs j viirdx — -Vcs- {2t 

'So "" 
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It remains to show that 

/ \WVi\'^Trdx= / [—^TT-Vi-— + VViTT]dS-VcS Vi1Tdx--VeS- (29) 

Jw Jso ^^ '^^ Jn' 2 

We note that L{vl) = 2viLvi + 2|Vt'ip = 2VesVi + 2\Vvi\'^. Thus 

/ IVvipvrdx = - / L{vl)7rdx — Vcs / Viirdx 
Jw 2 Jq, Jq, 

We took into account here that the integrand of the first term on the right vanishes at SK 
The last inequahty imphes fl29|) if the integral over 5*1 in the equality above is rewritten in 
terms of the integral over 5*0 using the first formula fl25l) and periodicity of the functions 
IT and Vi — s. The proof of the theorem is complete. Now let us justify the remark. 
From the divergence theorem and periodicity of tpi it follows that 

— — dx = / ipiUidS = / t/jiUidS. 
Since ni = ^ and the latter function on 5' is equal to — -^ (due to (l22l) ). we have 
—dx = / ipiTiidS = — / ipi^ — dS = — / ipi^ — dS. 



n' '^•s JdQ' Js^ ^^ J an' on 

Now we recall that Aipi = A(fi — s) = 0. Hence the latter formula together with the 
Green formula imply 

^'^idx = - [ \V^i\^dx. (30) 



'n' 9s jQi 

Further, if l^ = 0, then vr = l/|fi'|. From i;i = ^^i + s, d^S]) and (ED]) it follows that 

a^\n'\= [ \Vvi\^dx= / (l + 2^ + |V^i|')rfx= [ {l-\Vi^i\^)dx 
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Theorem 2.6. If V = then Ks = 0, and Kff = o-^V + 0{V^) as V ^ where 
(Tq = cr^|y=o > is the effective diffusivity atV = Q. In particular, if\V\ is small enough, 
then VcE 7^ when V ^ and Vcff is a monotone function ofV. 

Proof. From i^^ it follows that vr is an analytic function of V. Then fl25|) implies that 
Vcs is analytic in V, and ( 12^ implies that Vi and ipi depend on V analytically. We extend 
71 = Tjyr + Viii + O(y^) in the Taylor series at l^ = 0. Then (^^ leads to the following 
problem for tti: 

A7ri=0, xeil; -g^\s' = T^Ty ^i\si = t^i\soAt^i)s\si = [T^ijslso- 



Since ni = ^, only the factor — 1/1^2' | in the boundary condition makes the problem for 
TTi different from the problem for ip^ = ipi\v=o- The latter problem can be obtained from 
(EID if we put there V = V^s = 0. Hence tti = -^, i.e., tt = ^ - ^ + 0{V^). We put 
this into the second of relations (1251) and obtain that 

v„ = 1.(1 + i(|j>) + o(v^) = v(i - iiM^) + oiv-,. 



The last relation follows from fl30|) . It remains only to use the Remark after Theorem! 

D 

3 Periodic tubes with nearly separated dead zones. 

Consider the problem ([1]) with V = (without drift) in a periodic domain Q = 0,{e) C M'^ 
introduced in the previous section, that has the form of an e-independent periodic tube 
Qq with a periodic system of cavities connected to the tube Qq by narrow channels (see 
Fig. 3). In this section we'll assume that d > 3. We do not impose restrictions on the 
shape of channels, except the periodicity condition, smoothness of the boundary dQ{e) 
and the assumption that the channel enters the cell of periodicity Q'q of the main tube in 
an e neighborhood of some point Xq G dQ^. For simplicity we will assume that there is 
only one cavity and one channel on each period, but practically no changes are needed to 
extend all the arguments to the case of several cavities per period or multiple channels 
connecting the cavity with the tube. 




Figure 3: Periodic tube with nearly separated cavities (dead zones). 

Denote by ^'{e) = QQ[JQ[{e) one cell of periodicity of the domain Q{e), where Qq 
is the cell of periodicity of the main tube Qq and ^[{e) is the union of the cavity and 
the channel on this period. We assume that the cell Qq is defined by the restriction 
< s = Xi < 1. We denote by 5*0, 5*1 the parts of the boundary of Qq that belong to the 
planes through the points s = and s = 1, respectively. Let 5' = 9f2Q\[S'o IJ "S"!] be the 
lateral part of the boundary dfl'^ of the cell of the main tube. 

In order to describe the effective diffusivity a^ for the problem in Q{e) and the effective 
diffusivity ctq for the problem in the main tube Qq (without the cavities and channels), 
we introduce the following auxiliary problem 

Am = 0, X & Qq] u = s + ip, where ip is 1-periodic in s; u^\s' = 0. (31) 
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We will call function u the harmonic coordinate (for the domain Qq). This function was 
introduced in the previous section where it was denoted by vi (see (!22|) ). 

Theorem 3.1. Let V = 0. Then the following formulas are valid for the effective diffu- 
sivity (Tq in the case of a periodic tube Qq and for the effective diffusivity o"q(£:) in the case 
of the tube Qq with cavities: 

r„, iVul'^dx 

l"ol 

Remark. When Qq is a cylinder, formula (133|1 under some assumptions on the connecting 
channels can be found in [1], [1]. 

Proof. In the case of domain Qq, we will use notations u,ip for functions Vi,ipi, defined 
in fl2Tl) . fl22l) (compare fl22|l and fl3T!) ). and we preserve the original notations vi,ipi when 
the problem in ^{e) is considered. 

We note that vr is a constant when ^ = 0. From the normalization condition [j^, ndx = 
1) it follows that vr = I/I^qI for the problem in ^o and vr = l/|r2Q(£:)| for the problem 
in no(^)- Hence, formula ( 132|) follows from (126|) . In fact, ( l26l) was derived without any 
specific assumptions on the periodic domain Q, and it is valid for both domains Qq and 
n{e). Thus 

We are going to compare the functions u and Vi and derive f l33l) from flMl) and fl32l) . 

Let Xo G S*' be the point on the lateral side of the cell Q'q where the channel enters 
the main tube. Harmonic coordinates u and vi (for tubes Qq and f2(e), respectively) 
are defined up to arbitrary additive constants. We fix these constants assuming that 
u{xo) = t;i(xo) = 0. 

We fix a function a = a(e,x) G C°°(Qq) such that a = when \x — xo\ < 2e, a = 1 
when \x — xo\ > Se and |^ = on SK We consider function w = a{x)u{x) and extend it 
by zero in the channel and the cavity. Then w G C°° {Q' (e)) . 

Let us estimate Vi — w. Obviously, this difference satisfies the following relations in 

A(vi-w) = f ■.= 2VaVu + uAa, xenie); ^\~ ^' = 0, x E S\e), (35) 

on 

and it satisfies the periodicity condition on 5*0, 5*1. Function u is smooth and e-independent. 
Function Va has order 0{e~^), and its support belongs to a ball of radius 3e. Thus 
||VaVu||L2(n'(£)) = 0{e^~), £ — )■ 0. A similar estimate is valid for uAa since Aa = 
0{e~'^) and \u\ < Ce on the support of Aa. Hence ||/||L2(n'(e)) = 0{e^), e — )■ 0. We 
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also take into account that the support of / belongs to Q'q. Thus from the Green formula 
it follows that 

l|VK-w)||i^(f^,(,)) < C\\f\\L^(^n'ie))\\vi-w\\L2i%) = 0(e'^)||t;i-w||i2(n;,), £ -> 0. (36) 

Since domain fig does not depend on e and {vi — w){xo) = 0, we have \\vi — w\\l^(^q') < 
C||V(t'i — w^) 11^2(0;,). Hence from ( !36l) it follows that 

\\Vivi-w)\\L,in'(e)) = 0{e'^). s ^ 0. (37) 

Thus one can make the following changes in formula f l34|) with the accuracy of 0{e'^~'^): 
replace Vi by au, replace integration over ^'{e) by the integration over Q'q, and then drop 
a. In other words, one can replace the numerator in f lM|) by the numerator from fl32l) 
plus 0{e'^~'^). Then it remains only to use fl5^ and express the latter numerator through 
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